MMR-and-RMR-workflow-Clutch48
Elliott Schmidt
2024-06-27
Load libraries (packages)
System1 - Dell
Importing data from firesting for resting
preexperiment_date <- "21 March 2023 09 42AM/All"
postexperiment_date <- "21 March 2023 02 52PM/All"
##--- last fish run in trial ---##
experiment_date <- "21 March 2023 11 08AM/Oxygen"
experiment_date2 <- "21 March 2023 11 08AM/All"
firesting <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_date,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19)
Cycle_1 <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_date2,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE)
Cycle_last <-read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_date2,"slopes/Cycle_21.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) System2 - Asus
Importing data from firesting for resting
preexperiment_date_asus <- "21 March 2023 10 30AM/All"
postexperiment_date_asus <- "21 March 2023 03 16PM/All"
##--- last fish run in trial ---##
experiment_date_asus <- "21 March 2023 12 02PM/Oxygen"
experiment_date2_asus <- "21 March 2023 12 02PM/All"
firesting_asus <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_date_asus,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19)
Cycle_1_asus <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_date2_asus,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE)
Cycle_last_asus <-read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_date2_asus,"slopes/Cycle_21.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) Chamber volumes
chamber1_dell = 0.04650+0.00022
chamber2_dell = 0.04593+0.00022
chamber3_dell = 0.04977+0.00022
chamber4_dell = 0.04860+0.00022
chamber1_asus = 0.04565
chamber2_asus = 0.04573+0.00385
chamber3_asus = 0.04551+0.00322
chamber4_asus = 0.04791+0.00277
Date_tested="2023-03-21"
Clutch = "48"
Male = "CSUD013"
Female = "CSUD017"
Population = "Sudbury reef"
Tank =300
salinity =36
Date_analysed = Sys.Date() Replicates
1
Enter specimen data
Replicate = 1
mass = 0.0004185
chamber = "ch4"
Swim = "good/good"
chamber_vol = chamber4_dell
system1 = "Dell"
Notes=""
##--- time of trail ---##
experiment_mmr_date <- "21 March 2023 10 43AM/Oxygen"
experiment_mmr_date2 <- "21 March 2023 10 43AM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date2,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",preexperiment_date,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.002372409post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",postexperiment_date,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.003613005Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 2.16
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2$TIME ==Cycle_1$Time[1], firesting$TIME)
Tstart.dTIME=as.numeric(firesting2[Tstart.row, "dTIME"])
Tend.row=which(firesting2$TIME ==tail(Cycle_last$Time, n=1), firesting$TIME)
Tend.dTIME=as.numeric(firesting2[Tend.row, "dTIME"])
apoly_insp <- firesting2 |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 12 13 14 17 18 19 20 22 23 24
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 2.10
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=15,
measure=255,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...## plot.calc_rate.int: Plotting first 20 selected reps only. To plot others modify 'pos' input.
## -----------------------------------------##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from selected replicates...
## To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 0 rate(s) removed, 21 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 15 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 2 1 172.5409 -0.01282017 0.987 NA 492 717 5831.66
## 2: 5 1 197.8680 -0.01348117 0.998 NA 1923 2152 7451.98
## 3: 8 1 224.7284 -0.01402145 0.998 NA 3379 3610 9071.41
## 4: 11 1 250.1317 -0.01422998 0.995 NA 4841 5071 10693.06
## 5: 14 1 268.1448 -0.01387429 0.998 NA 6305 6536 12312.13
## 6: 21 1 334.8123 -0.01475019 0.994 NA 9723 9953 16093.19
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 6087.42 97.579 94.321 -0.01282017 -0.002537947 -0.01028222 -0.01028222
## 2: 7707.54 97.380 93.946 -0.01348117 -0.002739374 -0.01074180 -0.01074180
## 3: 9327.18 97.466 93.979 -0.01402145 -0.002940715 -0.01108073 -0.01108073
## 4: 10947.43 97.794 94.317 -0.01422998 -0.003142233 -0.01108774 -0.01108774
## 5: 12568.09 97.317 93.827 -0.01387429 -0.003343615 -0.01053067 -0.01053067
## 6: 16348.25 97.361 93.529 -0.01475019 -0.003813623 -0.01093657 -0.01093657
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04882 0.0004185 NA 36 27 1.013253 -0.1173138
## 2: %Air sec 0.04882 0.0004185 NA 36 27 1.013253 -0.1225573
## 3: %Air sec 0.04882 0.0004185 NA 36 27 1.013253 -0.1264243
## 4: %Air sec 0.04882 0.0004185 NA 36 27 1.013253 -0.1265043
## 5: %Air sec 0.04882 0.0004185 NA 36 27 1.013253 -0.1201485
## 6: %Air sec 0.04882 0.0004185 NA 36 27 1.013253 -0.1247796
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -280.3198 NA mgO2/hr/kg -280.3198
## 2: -292.8490 NA mgO2/hr/kg -292.8490
## 3: -302.0892 NA mgO2/hr/kg -302.0892
## 4: -302.2803 NA mgO2/hr/kg -302.2803
## 5: -287.0932 NA mgO2/hr/kg -287.0932
## 6: -298.1590 NA mgO2/hr/kg -298.1590
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 1 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.0004185 | ch4 | Dell | 0.04882 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 296.4942 | 0.1240828 | 0.9966 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 2.16
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 4 5 6 7 9 10 11 13 14 15 16 18 19 20 22 24 25 26
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 1.62
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 227.1635 -0.03354433 0.9681517 NA 21 75 3829.54
## 2: NA 2 227.1404 -0.03354149 0.9682797 NA 18 72 3826.19
## 3: NA 3 227.1085 -0.03353379 0.9683022 NA 17 71 3825.09
## 4: NA 4 227.0959 -0.03353114 0.9683464 NA 16 70 3824.00
## 5: NA 5 227.0906 -0.03352437 0.9679067 NA 22 76 3830.64
## ---
## 207: NA 207 153.0999 -0.01449044 0.9812756 NA 83 137 3899.61
## 208: NA 208 152.7321 -0.01439699 0.9815372 NA 84 138 3900.70
## 209: NA 209 152.5365 -0.01434759 0.9812678 NA 87 141 3904.18
## 210: NA 210 152.2454 -0.01427342 0.9814044 NA 86 140 3902.98
## 211: NA 211 152.2298 -0.01426927 0.9813829 NA 85 139 3901.86
## endtime oxy endoxy rate
## 1: 3889.54 98.580 96.806 -0.03354433
## 2: 3886.19 98.751 96.847 -0.03354149
## 3: 3885.09 98.780 96.890 -0.03353379
## 4: 3884.00 98.777 96.953 -0.03353114
## 5: 3890.64 98.568 96.778 -0.03352437
## ---
## 207: 3959.61 96.680 95.705 -0.01449044
## 208: 3960.70 96.646 95.718 -0.01439699
## 209: 3964.18 96.559 95.626 -0.01434759
## 210: 3962.98 96.596 95.660 -0.01427342
## 211: 3961.86 96.619 95.726 -0.01426927
##
## Regressions : 211 | Results : 211 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 211 adjusted rate(s):
## Rate : -0.03354433
## Adjustment : -0.002372409
## Adjusted Rate : -0.03117192
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 19 rate(s) removed, 192 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 191 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 227.1635 -0.03354433 0.9681517 NA 21 75 3829.54
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 3889.54 98.58 96.806 -0.03354433 -0.002372409 -0.03117192 -0.03117192
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04882 0.0004185 NA 36 27 1.013253 -0.3556525
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -849.8267 NA mgO2/hr/kg -849.8267
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 1 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.0004185 | ch4 | Dell | 0.04882 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 296.4942 | 0.1240828 | 0.9966 | 849.8267 | 0.3556525 | 0.9681517 | 553.3325 | 0.2315697 |
Exporting data
## Rows: 327 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.2
Enter specimen data
Replicate = 2
mass = 0.0004897
chamber = "ch3"
Swim = "good/good"
chamber_vol = chamber3_dell
system1 = "Dell"
Notes="check max"
##--- time of trail ---##
experiment_mmr_date <- "21 March 2023 10 55AM/Oxygen"
experiment_mmr_date2 <- "21 March 2023 10 55AM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date2,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",preexperiment_date,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.001644409post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",postexperiment_date,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.002810159Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 2.16
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2$TIME ==Cycle_1$Time[1], firesting$TIME)
Tstart.dTIME=as.numeric(firesting2[Tstart.row, "dTIME"])
Tend.row=which(firesting2$TIME ==tail(Cycle_last$Time, n=1), firesting$TIME)
Tend.dTIME=as.numeric(firesting2[Tend.row, "dTIME"])
apoly_insp <- firesting2 |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 12 13 14 17 18 19 20 22 23 24
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 2.10
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=15,
measure=255,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...## plot.calc_rate.int: Plotting first 20 selected reps only. To plot others modify 'pos' input.
## -----------------------------------------##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from selected replicates...
## To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 2 rate(s) removed, 19 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 13 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 1 1 142.7094 -0.008412286 0.957 NA 15 241 5292.92
## 2: 2 1 162.3973 -0.011390172 0.988 NA 492 717 5831.66
## 3: 10 1 197.6579 -0.009838487 0.987 NA 4353 4583 10151.81
## 4: 13 1 243.0901 -0.012325593 0.984 NA 5816 6047 11771.57
## 5: 14 1 239.0211 -0.011453449 0.985 NA 6305 6536 12312.13
## 6: 21 1 295.0422 -0.012252267 0.989 NA 9723 9953 16093.19
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 5547.37 98.227 96.203 -0.008412286 -0.001736947 -0.006675338 -0.006675338
## 2: 6087.42 95.778 92.949 -0.011390172 -0.001799960 -0.009590212 -0.009590212
## 3: 10407.22 97.739 95.266 -0.009838487 -0.002304621 -0.007533866 -0.007533866
## 4: 12026.84 97.747 94.704 -0.012325593 -0.002493834 -0.009831759 -0.009831759
## 5: 12568.09 97.849 94.736 -0.011453449 -0.002557023 -0.008896426 -0.008896426
## 6: 16348.25 98.012 94.756 -0.012252267 -0.002998675 -0.009253592 -0.009253592
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04999 0.0004897 NA 36 27 1.013253 -0.07798676
## 2: %Air sec 0.04999 0.0004897 NA 36 27 1.013253 -0.11204070
## 3: %Air sec 0.04999 0.0004897 NA 36 27 1.013253 -0.08801679
## 4: %Air sec 0.04999 0.0004897 NA 36 27 1.013253 -0.11486265
## 5: %Air sec 0.04999 0.0004897 NA 36 27 1.013253 -0.10393533
## 6: %Air sec 0.04999 0.0004897 NA 36 27 1.013253 -0.10810803
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -159.2542 NA mgO2/hr/kg -159.2542
## 2: -228.7946 NA mgO2/hr/kg -228.7946
## 3: -179.7361 NA mgO2/hr/kg -179.7361
## 4: -234.5572 NA mgO2/hr/kg -234.5572
## 5: -212.2429 NA mgO2/hr/kg -212.2429
## 6: -220.7638 NA mgO2/hr/kg -220.7638
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 2 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.0004897 | ch3 | Dell | 0.04999 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 215.2189 | 0.1053927 | 0.9866 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 2.16
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 15 17 18 19 20 21 24
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 1.55
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 236.7964 -0.03019152 0.9883818 NA 216 271 4747.34
## 2: NA 2 236.4521 -0.03011963 0.9882053 NA 215 270 4746.24
## 3: NA 3 236.2043 -0.03006780 0.9881396 NA 214 269 4745.16
## 4: NA 4 235.8133 -0.02998642 0.9876818 NA 213 268 4744.06
## 5: NA 5 234.9984 -0.02981638 0.9864553 NA 212 267 4742.95
## ---
## 212: NA 212 159.2162 -0.01345080 0.9196025 NA 5 60 4512.50
## 213: NA 213 157.1323 -0.01299176 0.9187848 NA 4 59 4511.40
## 214: NA 214 155.2579 -0.01257912 0.9114893 NA 3 58 4510.30
## 215: NA 215 155.0997 -0.01254432 0.9109537 NA 2 57 4509.21
## 216: NA 216 154.7489 -0.01246618 0.9128784 NA 1 56 4508.11
## endtime oxy endoxy rate
## 1: 4807.34 93.438 91.660 -0.03019152
## 2: 4806.24 93.491 91.683 -0.03011963
## 3: 4805.16 93.453 91.769 -0.03006780
## 4: 4804.06 93.434 91.799 -0.02998642
## 5: 4802.95 93.454 91.750 -0.02981638
## ---
## 212: 4572.50 98.427 97.562 -0.01345080
## 213: 4571.40 98.395 97.649 -0.01299176
## 214: 4570.30 98.516 97.766 -0.01257912
## 215: 4569.21 98.633 97.656 -0.01254432
## 216: 4568.11 98.691 97.663 -0.01246618
##
## Regressions : 216 | Results : 216 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 216 adjusted rate(s):
## Rate : -0.03019152
## Adjustment : -0.001644409
## Adjusted Rate : -0.02854711
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 11 rate(s) removed, 205 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 204 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 236.7964 -0.03019152 0.9883818 NA 216 271 4747.34
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 4807.34 93.438 91.66 -0.03019152 -0.001644409 -0.02854711 -0.02854711
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04999 0.0004897 NA 36 27 1.013253 -0.3335107
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -681.0512 NA mgO2/hr/kg -681.0512
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 2 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.0004897 | ch3 | Dell | 0.04999 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 215.2189 | 0.1053927 | 0.9866 | 681.0512 | 0.3335107 | 0.9883818 | 465.8322 | 0.2281181 | check max |
Exporting data
## Rows: 328 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.3
Enter specimen data
Replicate = 3
mass = 0.0004472
chamber = "ch2"
Swim = "good/good"
chamber_vol = chamber2_dell
system1 = "Dell"
Notes="check max"
##--- time of trail ---##
experiment_mmr_date <- "21 March 2023 11 08AM/Oxygen"
experiment_mmr_date2 <- "21 March 2023 11 08AM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",experiment_mmr_date2,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",preexperiment_date,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.001547705post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Dell/Experiment_",postexperiment_date,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.002849508Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 2.16
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2$TIME ==Cycle_1$Time[1], firesting$TIME)
Tstart.dTIME=as.numeric(firesting2[Tstart.row, "dTIME"])
Tend.row=which(firesting2$TIME ==tail(Cycle_last$Time, n=1), firesting$TIME)
Tend.dTIME=as.numeric(firesting2[Tend.row, "dTIME"])
apoly_insp <- firesting2 |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 12 13 14 17 18 19 20 22 23 24
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 2.10
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=15,
measure=255,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...## plot.calc_rate.int: Plotting first 20 selected reps only. To plot others modify 'pos' input.
## -----------------------------------------##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from selected replicates...
## To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 2 rate(s) removed, 19 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 13 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 2 1 165.3543 -0.011669067 0.987 NA 492 717 5831.66
## 2: 6 1 174.2539 -0.009629711 0.979 NA 2407 2637 7992.39
## 3: 10 1 239.0554 -0.013997436 0.984 NA 4353 4583 10151.81
## 4: 11 1 245.8077 -0.013860354 0.964 NA 4841 5071 10693.06
## 5: 12 1 251.8677 -0.013741514 0.989 NA 5329 5560 11232.98
## 6: 13 1 263.0261 -0.014088380 0.979 NA 5816 6047 11771.57
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 6087.42 97.456 94.383 -0.011669067 -0.001721410 -0.009947657 -0.009947657
## 2: 8247.88 97.233 94.741 -0.009629711 -0.002003269 -0.007626442 -0.007626442
## 3: 10407.22 96.855 93.265 -0.013997436 -0.002284969 -0.011712467 -0.011712467
## 4: 10947.43 97.302 93.741 -0.013860354 -0.002355510 -0.011504844 -0.011504844
## 5: 11488.04 97.439 93.832 -0.013741514 -0.002425989 -0.011315524 -0.011315524
## 6: 12026.84 96.911 93.409 -0.014088380 -0.002496264 -0.011592116 -0.011592116
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04615 0.0004472 NA 36 27 1.013253 -0.10728944
## 2: %Air sec 0.04615 0.0004472 NA 36 27 1.013253 -0.08225422
## 3: %Air sec 0.04615 0.0004472 NA 36 27 1.013253 -0.12632362
## 4: %Air sec 0.04615 0.0004472 NA 36 27 1.013253 -0.12408433
## 5: %Air sec 0.04615 0.0004472 NA 36 27 1.013253 -0.12204244
## 6: %Air sec 0.04615 0.0004472 NA 36 27 1.013253 -0.12502559
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -239.9138 NA mgO2/hr/kg -239.9138
## 2: -183.9316 NA mgO2/hr/kg -183.9316
## 3: -282.4768 NA mgO2/hr/kg -282.4768
## 4: -277.4694 NA mgO2/hr/kg -277.4694
## 5: -272.9035 NA mgO2/hr/kg -272.9035
## 6: -279.5742 NA mgO2/hr/kg -279.5742
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 3 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.0004472 | ch2 | Dell | 0.04615 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 270.4675 | 0.1209531 | 0.9806 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 2.16
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 3 4 5 6 7 8 9 10 11 12 13 15 16 17 18 19 20 21 22
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.08 1.35
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 312.5412 -0.04004428 0.9932208 NA 165 220 5463.45
## 2: NA 2 312.3193 -0.04000417 0.9931421 NA 164 219 5462.36
## 3: NA 3 311.7062 -0.03989301 0.9927997 NA 163 218 5461.22
## 4: NA 4 311.3754 -0.03983173 0.9926901 NA 166 221 5464.64
## 5: NA 5 310.7154 -0.03971302 0.9922517 NA 162 217 5460.11
## ---
## 208: NA 208 205.0894 -0.02013935 0.9608066 NA 5 60 5282.49
## 209: NA 209 201.8465 -0.01952945 0.9506266 NA 4 59 5281.38
## 210: NA 210 198.8320 -0.01896245 0.9398451 NA 3 58 5280.27
## 211: NA 211 196.6297 -0.01854791 0.9352285 NA 2 57 5279.18
## 212: NA 212 194.4936 -0.01814555 0.9327046 NA 1 56 5278.05
## endtime oxy endoxy rate
## 1: 5523.45 93.736 91.391 -0.04004428
## 2: 5522.36 93.793 91.440 -0.04000417
## 3: 5521.22 93.890 91.467 -0.03989301
## 4: 5524.64 93.713 91.424 -0.03983173
## 5: 5520.11 93.905 91.449 -0.03971302
## ---
## 208: 5342.49 98.455 97.415 -0.02013935
## 209: 5341.38 98.463 97.468 -0.01952945
## 210: 5340.27 98.565 97.488 -0.01896245
## 211: 5339.18 98.608 97.491 -0.01854791
## 212: 5338.05 98.558 97.535 -0.01814555
##
## Regressions : 212 | Results : 212 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 212 adjusted rate(s):
## Rate : -0.04004428
## Adjustment : -0.001547705
## Adjusted Rate : -0.03849658
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 21 rate(s) removed, 191 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 190 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 312.5412 -0.04004428 0.9932208 NA 165 220 5463.45
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 5523.45 93.736 91.391 -0.04004428 -0.001547705 -0.03849658 -0.03849658
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04615 0.0004472 NA 36 27 1.013253 -0.4152009
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -928.4456 NA mgO2/hr/kg -928.4456
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 3 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.0004472 | ch2 | Dell | 0.04615 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 270.4675 | 0.1209531 | 0.9806 | 928.4456 | 0.4152009 | 0.9932208 | 657.9781 | 0.2942478 | check max |
Exporting data
## Rows: 329 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.5
Enter specimen data
Replicate = 5
mass = 0.0004890
chamber = "ch4"
Swim = "good/good"
chamber_vol = chamber4_asus
system1 = "Asus"
Notes=""
##--- time of trail ---##
experiment_mmr_date_asus <- "21 March 2023 11 24AM/Oxygen"
experiment_mmr_date2_asus <- "21 March 2023 11 24AM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date_asus,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date2_asus,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",preexperiment_date_asus,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.0005512114post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",postexperiment_date_asus,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.002265624Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 4.70
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2_asus$TIME ==Cycle_1_asus$Time[1], firesting_asus$TIME)
Tstart.dTIME=as.numeric(firesting2_asus[Tstart.row, "dTIME"])
Tend.row=which(firesting2_asus$TIME ==tail(Cycle_last_asus$Time, n=1), firesting_asus$TIME)
Tend.dTIME=as.numeric(firesting2_asus[Tend.row, "dTIME"])
apoly_insp <- firesting2_asus |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 6 7 8 9 10 11 12 13 15 16 17 18 19 21 22 23
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 4.70
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=15,
measure=255,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...## plot.calc_rate.int: Plotting first 20 selected reps only. To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2_asus$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 0 rate(s) removed, 21 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 15 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 5 1 204.7827 -0.01355224 0.984 NA 1610 1799 7809.60
## 2: 17 1 322.5826 -0.01566114 0.995 NA 6412 6599 14290.21
## 3: 18 1 338.0205 -0.01610380 0.994 NA 6810 6959 14829.55
## 4: 19 1 339.6235 -0.01564471 0.982 NA 7129 7318 15369.25
## 5: 20 1 357.0156 -0.01621662 0.992 NA 7530 7719 15910.30
## 6: 21 1 377.2550 -0.01689898 0.991 NA 7928 8115 16450.41
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 8064.94 98.608 95.439 -0.01355224 -0.001371708 -0.01218054 -0.01218054
## 2: 14544.97 98.595 94.669 -0.01566114 -0.002737069 -0.01292407 -0.01292407
## 3: 15084.89 98.982 95.048 -0.01610380 -0.002850765 -0.01325303 -0.01325303
## 4: 15624.45 98.838 94.979 -0.01564471 -0.002964462 -0.01268025 -0.01268025
## 5: 16165.84 98.568 94.824 -0.01621662 -0.003078493 -0.01313813 -0.01313813
## 6: 16705.01 98.970 94.858 -0.01689898 -0.003192192 -0.01370679 -0.01370679
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.05068 0.000489 NA 36 27 1.013253 -0.1442672
## 2: %Air sec 0.05068 0.000489 NA 36 27 1.013253 -0.1530737
## 3: %Air sec 0.05068 0.000489 NA 36 27 1.013253 -0.1569699
## 4: %Air sec 0.05068 0.000489 NA 36 27 1.013253 -0.1501858
## 5: %Air sec 0.05068 0.000489 NA 36 27 1.013253 -0.1556090
## 6: %Air sec 0.05068 0.000489 NA 36 27 1.013253 -0.1623442
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -295.0249 NA mgO2/hr/kg -295.0249
## 2: -313.0341 NA mgO2/hr/kg -313.0341
## 3: -321.0018 NA mgO2/hr/kg -321.0018
## 4: -307.1284 NA mgO2/hr/kg -307.1284
## 5: -318.2187 NA mgO2/hr/kg -318.2187
## 6: -331.9922 NA mgO2/hr/kg -331.9922
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 5 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.000489 | ch4 | Asus | 0.05068 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 318.2751 | 0.1556365 | 0.9908 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 3.85
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch4
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 6 7 8 9 10 11 12 13 14 15 17 18 19 20 22 23
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 1.67
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 236.3819 -0.04015514 0.9883623 NA 71 116 3466.29
## 2: NA 2 236.2685 -0.04012445 0.9881641 NA 70 115 3464.95
## 3: NA 3 235.9882 -0.04004574 0.9877304 NA 69 114 3463.62
## 4: NA 4 235.5708 -0.03992168 0.9871064 NA 72 117 3467.62
## 5: NA 5 235.1124 -0.03979537 0.9872920 NA 68 113 3462.26
## ---
## 174: NA 174 161.1504 -0.01842480 0.9882756 NA 1 46 3371.55
## 175: NA 175 161.1409 -0.01842157 0.9885729 NA 8 53 3381.00
## 176: NA 176 161.0735 -0.01840171 0.9886381 NA 3 48 3374.25
## 177: NA 177 161.0205 -0.01838660 0.9887233 NA 5 50 3376.96
## 178: NA 178 160.9550 -0.01836719 0.9887470 NA 4 49 3375.62
## endtime oxy endoxy rate
## 1: 3526.29 97.085 94.911 -0.04015514
## 2: 3524.95 97.139 94.931 -0.04012445
## 3: 3523.62 97.237 94.855 -0.04004574
## 4: 3527.62 97.113 94.910 -0.03992168
## 5: 3522.26 97.253 94.899 -0.03979537
## ---
## 174: 3431.55 99.040 97.908 -0.01842480
## 175: 3441.00 98.889 97.765 -0.01842157
## 176: 3434.25 99.032 97.937 -0.01840171
## 177: 3436.96 98.953 97.810 -0.01838660
## 178: 3435.62 99.002 97.837 -0.01836719
##
## Regressions : 178 | Results : 178 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 178 adjusted rate(s):
## Rate : -0.04015514
## Adjustment : -0.0005512114
## Adjusted Rate : -0.03960393
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 0 rate(s) removed, 178 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 177 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 236.3819 -0.04015514 0.9883623 NA 71 116 3466.29
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 3526.29 97.085 94.911 -0.04015514 -0.0005512114 -0.03960393 -0.03960393
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.05068 0.000489 NA 36 27 1.013253 -0.4690719
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -959.2472 NA mgO2/hr/kg -959.2472
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 5 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.000489 | ch4 | Asus | 0.05068 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 318.2751 | 0.1556365 | 0.9908 | 959.2472 | 0.4690719 | 0.9883623 | 640.9721 | 0.3134354 | ||
| ### Expor | ting data |
## Rows: 330 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.6
Enter specimen data
Replicate = 6
mass = 0.0003880
chamber = "ch3"
Swim = "good/good"
chamber_vol = chamber3_asus
system1 = "Asus"
Notes=""
##--- time of trail ---##
experiment_mmr_date_asus <- "21 March 2023 11 38AM/Oxygen"
experiment_mmr_date2_asus <- "21 March 2023 11 38AM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date_asus,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date2_asus,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",preexperiment_date_asus,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] 0.0001501796post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",postexperiment_date_asus,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.001770659Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 4.70
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2_asus$TIME ==Cycle_1_asus$Time[1], firesting_asus$TIME)
Tstart.dTIME=as.numeric(firesting2_asus[Tstart.row, "dTIME"])
Tend.row=which(firesting2_asus$TIME ==tail(Cycle_last_asus$Time, n=1), firesting_asus$TIME)
Tend.dTIME=as.numeric(firesting2_asus[Tend.row, "dTIME"])
apoly_insp <- firesting2_asus |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 6 7 8 9 10 11 12 13 15 16 17 18 19 21 22 23
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 4.70
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=15,
measure=245,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...## plot.calc_rate.int: Plotting first 20 selected reps only. To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## Warning: adjust_rate: background rates in 'by' and 'by2' differ in sign (i.e. one is +ve, one is -ve).
## Ensure this is correct. The 'linear' adjustment has been performed regardless.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 0 rate(s) removed, 21 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 15 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 2 1 163.1779 -0.01046895 0.991 NA 411 592 6189.38
## 2: 3 1 173.0462 -0.01105577 0.989 NA 811 993 6730.21
## 3: 14 1 249.3413 -0.01195703 0.993 NA 5211 5393 12669.04
## 4: 15 1 261.8189 -0.01238582 0.998 NA 5612 5794 13209.79
## 5: 18 1 287.3304 -0.01272858 0.995 NA 6810 6954 14829.55
## 6: 19 1 287.9160 -0.01232833 0.998 NA 7129 7311 15369.25
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 6434.11 98.373 95.669 -0.01046895 -0.0003853853 -0.010083568 -0.010083568
## 2: 6975.94 98.452 95.919 -0.01105577 -0.0005131729 -0.010542594 -0.010542594
## 3: 12914.46 97.964 95.018 -0.01195703 -0.0019150696 -0.010041961 -0.010041961
## 4: 13455.42 98.248 95.220 -0.01238582 -0.0020427450 -0.010343079 -0.010343079
## 5: 15074.95 98.590 95.311 -0.01272858 -0.0024250819 -0.010303496 -0.010303496
## 6: 15615.01 98.365 95.430 -0.01232833 -0.0025525272 -0.009775798 -0.009775798
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04873 0.000388 NA 36 27 1.013253 -0.1148352
## 2: %Air sec 0.04873 0.000388 NA 36 27 1.013253 -0.1200628
## 3: %Air sec 0.04873 0.000388 NA 36 27 1.013253 -0.1143614
## 4: %Air sec 0.04873 0.000388 NA 36 27 1.013253 -0.1177906
## 5: %Air sec 0.04873 0.000388 NA 36 27 1.013253 -0.1173398
## 6: %Air sec 0.04873 0.000388 NA 36 27 1.013253 -0.1113302
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -295.9671 NA mgO2/hr/kg -295.9671
## 2: -309.4402 NA mgO2/hr/kg -309.4402
## 3: -294.7458 NA mgO2/hr/kg -294.7458
## 4: -303.5841 NA mgO2/hr/kg -303.5841
## 5: -302.4223 NA mgO2/hr/kg -302.4223
## 6: -286.9336 NA mgO2/hr/kg -286.9336
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 6 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.000388 | ch3 | Asus | 0.04873 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 301.2319 | 0.116878 | 0.9932 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 3.85
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch3
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 2 3 4 5 6 7 8 9 10 12 13 14 15 16 17 18 19 22 23 24
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.33 1.45
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 237.5387 -0.03276340 0.9754182 NA 52 97 4264.53
## 2: NA 2 237.3124 -0.03271183 0.9750119 NA 51 96 4263.20
## 3: NA 3 237.2240 -0.03268862 0.9746650 NA 53 98 4265.89
## 4: NA 4 237.0460 -0.03265098 0.9745042 NA 50 95 4261.86
## 5: NA 5 236.7611 -0.03257929 0.9735237 NA 54 99 4267.34
## ---
## 174: NA 174 148.6148 -0.01231663 0.9661869 NA 177 222 4433.28
## 175: NA 175 148.6042 -0.01231457 0.9658312 NA 175 220 4430.62
## 176: NA 176 148.2341 -0.01223140 0.9647576 NA 173 218 4427.90
## 177: NA 177 148.1148 -0.01220491 0.9649038 NA 174 219 4429.26
## 178: NA 178 147.4355 -0.01205263 0.9666000 NA 178 223 4434.65
## endtime oxy endoxy rate
## 1: 4324.53 97.713 95.949 -0.03276340
## 2: 4323.20 97.742 95.993 -0.03271183
## 3: 4325.89 97.666 95.976 -0.03268862
## 4: 4321.86 97.786 95.957 -0.03265098
## 5: 4327.34 97.625 95.954 -0.03257929
## ---
## 174: 4493.28 94.035 93.311 -0.01231663
## 175: 4490.62 94.018 93.288 -0.01231457
## 176: 4487.90 94.175 93.325 -0.01223140
## 177: 4489.26 94.099 93.303 -0.01220491
## 178: 4494.65 94.076 93.303 -0.01205263
##
## Regressions : 178 | Results : 178 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 178 adjusted rate(s):
## Rate : -0.0327634
## Adjustment : 0.0001501796
## Adjusted Rate : -0.03291358
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 0 rate(s) removed, 178 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 177 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 237.5387 -0.0327634 0.9754182 NA 52 97 4264.53
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 4324.53 97.713 95.949 -0.0327634 0.0001501796 -0.03291358 -0.03291358
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04873 0.000388 NA 36 27 1.013253 -0.3748314
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -966.0604 NA mgO2/hr/kg -966.0604
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 6 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.000388 | ch3 | Asus | 0.04873 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 301.2319 | 0.116878 | 0.9932 | 966.0604 | 0.3748314 | 0.9754182 | 664.8285 | 0.2579535 | ||
| ### Expor | ting data |
## Rows: 331 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.7
Enter specimen data
Replicate = 7
mass = 0.0004910
chamber = "ch2"
Swim = "good/good"
chamber_vol = chamber2_asus
system1 = "Asus"
Notes="check max"
##--- time of trail ---##
experiment_mmr_date_asus <- "21 March 2023 11 50AM/Oxygen"
experiment_mmr_date2_asus <- "21 March 2023 11 50AM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date_asus,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date2_asus,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",preexperiment_date_asus,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] 0.0005094143post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",postexperiment_date_asus,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.002201281Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 4.70
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2_asus$TIME ==Cycle_1_asus$Time[1], firesting_asus$TIME)
Tstart.dTIME=as.numeric(firesting2_asus[Tstart.row, "dTIME"])
Tend.row=which(firesting2_asus$TIME ==tail(Cycle_last_asus$Time, n=1), firesting_asus$TIME)
Tend.dTIME=as.numeric(firesting2_asus[Tend.row, "dTIME"])
apoly_insp <- firesting2_asus |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 6 7 8 9 10 11 12 13 15 16 17 18 19 21 22 23
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 4.70
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=15,
measure=245,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...## plot.calc_rate.int: Plotting first 20 selected reps only. To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## Warning: adjust_rate: background rates in 'by' and 'by2' differ in sign (i.e. one is +ve, one is -ve).
## Ensure this is correct. The 'linear' adjustment has been performed regardless.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 0 rate(s) removed, 21 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 15 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 12 1 264.4144 -0.01435843 0.994 NA 4410 4592 11589.17
## 2: 17 1 315.7471 -0.01523008 0.995 NA 6412 6592 14290.21
## 3: 18 1 305.1368 -0.01394527 0.997 NA 6810 6954 14829.55
## 4: 19 1 325.9144 -0.01480435 0.990 NA 7129 7311 15369.25
## 5: 20 1 328.4002 -0.01445795 0.991 NA 7530 7711 15910.30
## 6: 21 1 349.9226 -0.01527885 0.995 NA 7928 8108 16450.41
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 11834.78 97.975 94.455 -0.01435843 -0.002045367 -0.01231306 -0.01231306
## 2: 14535.37 97.927 94.207 -0.01523008 -0.002945095 -0.01228499 -0.01228499
## 3: 15074.95 98.232 94.878 -0.01394527 -0.003124806 -0.01082046 -0.01082046
## 4: 15615.01 98.124 94.608 -0.01480435 -0.003304657 -0.01149969 -0.01149969
## 5: 16154.94 98.127 94.616 -0.01445795 -0.003484712 -0.01097324 -0.01097324
## 6: 16695.55 98.401 94.714 -0.01527885 -0.003664723 -0.01161413 -0.01161413
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04958 0.000491 NA 36 27 1.013253 -0.1426714
## 2: %Air sec 0.04958 0.000491 NA 36 27 1.013253 -0.1423462
## 3: %Air sec 0.04958 0.000491 NA 36 27 1.013253 -0.1253767
## 4: %Air sec 0.04958 0.000491 NA 36 27 1.013253 -0.1332469
## 5: %Air sec 0.04958 0.000491 NA 36 27 1.013253 -0.1271469
## 6: %Air sec 0.04958 0.000491 NA 36 27 1.013253 -0.1345729
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -290.5732 NA mgO2/hr/kg -290.5732
## 2: -289.9107 NA mgO2/hr/kg -289.9107
## 3: -255.3496 NA mgO2/hr/kg -255.3496
## 4: -271.3786 NA mgO2/hr/kg -271.3786
## 5: -258.9550 NA mgO2/hr/kg -258.9550
## 6: -274.0793 NA mgO2/hr/kg -274.0793
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 7 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.000491 | ch2 | Asus | 0.04958 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 276.9794 | 0.1359969 | 0.993 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 3.85
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch2
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 2 3 4 5 6 7 8 9 10 13 14 15 16 17 18 19 20 23 24 25
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 1.49
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 254.8280 -0.03143283 0.9882998 NA 133 178 5075.03
## 2: NA 2 254.8035 -0.03142877 0.9882770 NA 132 177 5073.68
## 3: NA 3 253.7045 -0.03121192 0.9866624 NA 134 179 5076.38
## 4: NA 4 253.4969 -0.03117325 0.9870818 NA 131 176 5072.35
## 5: NA 5 252.3308 -0.03094207 0.9846076 NA 135 180 5077.73
## ---
## 174: NA 174 154.3552 -0.01182779 0.8334557 NA 163 208 5115.54
## 175: NA 175 155.3545 -0.01166520 0.9190509 NA 4 49 4900.79
## 176: NA 176 153.5414 -0.01129655 0.9306124 NA 3 48 4899.42
## 177: NA 177 151.8526 -0.01095305 0.9425320 NA 2 47 4898.09
## 178: NA 178 151.0758 -0.01079437 0.9501698 NA 1 46 4896.76
## endtime oxy endoxy rate
## 1: 5135.03 95.241 93.501 -0.03143283
## 2: 5133.68 95.243 93.462 -0.03142877
## 3: 5136.38 95.234 93.534 -0.03121192
## 4: 5132.35 95.231 93.487 -0.03117325
## 5: 5137.73 95.200 93.515 -0.03094207
## ---
## 174: 5175.54 94.091 92.928 -0.01182779
## 175: 4960.79 98.209 97.308 -0.01166520
## 176: 4959.42 98.219 97.349 -0.01129655
## 177: 4958.09 98.291 97.400 -0.01095305
## 178: 4956.76 98.308 97.465 -0.01079437
##
## Regressions : 178 | Results : 178 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 178 adjusted rate(s):
## Rate : -0.03143283
## Adjustment : 0.0005094143
## Adjusted Rate : -0.03194224
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 44 rate(s) removed, 134 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 133 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 254.828 -0.03143283 0.9882998 NA 133 178 5075.03
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 5135.03 95.241 93.501 -0.03143283 0.0005094143 -0.03194224 -0.03194224
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04958 0.000491 NA 36 27 1.013253 -0.3701148
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -753.7979 NA mgO2/hr/kg -753.7979
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 7 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.000491 | ch2 | Asus | 0.04958 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 276.9794 | 0.1359969 | 0.993 | 753.7979 | 0.3701148 | 0.9882998 | 476.8185 | 0.2341179 | check max | |
| ### Expor | ting data |
## Rows: 332 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.8
Enter specimen data
Replicate = 8
mass = 0.0005137
chamber = "ch1"
Swim = "good/good"
chamber_vol = chamber1_asus
system1 = "Asus"
Notes=""
##--- time of trail ---##
experiment_mmr_date_asus <- "21 March 2023 12 02PM/Oxygen"
experiment_mmr_date2_asus <- "21 March 2023 12 02PM/All"
firesting_mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date_asus,"data raw/Firesting.txt"),
delim = "\t", escape_double = FALSE,
col_types = cols(`Time (HH:MM:SS)` = col_time(format = "%H:%M:%S"),
`Time (s)` = col_number(), Ch1...5 = col_number(),
Ch2...6 = col_number(), Ch3...7 = col_number(),
Ch4...8 = col_number()), trim_ws = TRUE,
skip = 19) ## New names:
## • `Ch1` -> `Ch1...5`
## • `Ch2` -> `Ch2...6`
## • `Ch3` -> `Ch3...7`
## • `Ch4` -> `Ch4...8`
## • `Ch 1` -> `Ch 1...9`
## • `Ch 2` -> `Ch 2...10`
## • `Ch 3` -> `Ch 3...11`
## • `Ch 4` -> `Ch 4...12`
## • `('C)` -> `('C)...15`
## • `('C)` -> `('C)...16`
## • `Ch 1` -> `Ch 1...18`
## • `Ch 2` -> `Ch 2...19`
## • `Ch 3` -> `Ch 3...20`
## • `Ch 4` -> `Ch 4...21`
## • `Ch1` -> `Ch1...22`
## • `Ch2` -> `Ch2...23`
## • `Ch3` -> `Ch3...24`
## • `Ch4` -> `Ch4...25`
## • `Ch1` -> `Ch1...26`
## • `Ch2` -> `Ch2...27`
## • `Ch3` -> `Ch3...28`
## • `Ch4` -> `Ch4...29`
## • `` -> `...31`## Warning: One or more parsing issues, call `problems()` on your data frame for details,
## e.g.:
## dat <- vroom(...)
## problems(dat)Cycle_1.mmr <- read_delim(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",experiment_mmr_date2_asus,"slopes/Cycle_1.txt"),
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
`Seconds from start for linreg` = col_number(),
`ch1 po2` = col_number(), `ch2 po2` = col_number(),
`ch3 po2` = col_number(), `ch4 po2` = col_number(),
...8 = col_skip()), trim_ws = TRUE) ## New names:
## • `` -> `...8`Background rates
Pre-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",preexperiment_date_asus,"slopes"))
pre_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
pre_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_pre1 <- pre_cycle1 %>% calc_rate.bg()##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] 0.0008872263post-experiment
setwd(paste("C:/Users/jc527762/OneDrive - James Cook University/PhD dissertation/Data/2023/Resp_backup/2023_Resp/Asus/Experiment_",postexperiment_date_asus,"slopes"))
post_cycle1 <- read_delim("./Cycle_1.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle2 <- read_delim("./Cycle_2.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
post_cycle3 <- read_delim("./Cycle_3.txt",
delim = ";", escape_double = FALSE, col_types = cols(Time = col_time(format = "%H:%M:%S"),
...8 = col_skip()), trim_ws = TRUE) %>%
rename(dTIME = `Seconds from start for linreg`,
ch1 =`ch1 po2`,
ch2 =`ch2 po2`,
ch3 =`ch3 po2`,
ch4 =`ch4 po2`) %>%
select(c("Time",chamber))
bg_post1 <- post_cycle1 %>% calc_rate.bg() ##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------##
## # plot.calc_rate.bg # -------------------
## plot.calc_rate.bg: Plotting all 1 background rates ...
## -----------------------------------------## [1] -0.00226215Resting metabolic rate
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 4.70
## -----------------------------------------
#### subset data
Tstart.row=which(firesting2_asus$TIME ==Cycle_1_asus$Time[1], firesting_asus$TIME)
Tstart.dTIME=as.numeric(firesting2_asus[Tstart.row, "dTIME"])
Tend.row=which(firesting2_asus$TIME ==tail(Cycle_last_asus$Time, n=1), firesting_asus$TIME)
Tend.dTIME=as.numeric(firesting2_asus[Tend.row, "dTIME"])
apoly_insp <- firesting2_asus |>
subset_data(from=Tstart.dTIME,
to=Tend.dTIME,
by="time")
apoly_insp <- inspect(apoly_insp, time=1, oxygen=2)##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 6 7 8 9 10 11 12 13 15 16 17 18 19 21 22 23
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 4.70
## -----------------------------------------
Extract rates
apoly_cr.int <- calc_rate.int(apoly_insp,
starts=(195+45+300),
wait=15,
measure=245,
by="time",
plot=TRUE) ##
## # plot.calc_rate.int # ------------------
## plot.calc_rate.int: Plotting rate from all replicates ...## plot.calc_rate.int: Plotting first 20 selected reps only. To plot others modify 'pos' input.
## -----------------------------------------adjust rates for background
apoly_cr.int_adj <- adjust_rate(apoly_cr.int,
by = bg_pre,
by2 = bg_post,
time_by = Tstart.row,
time_by2 = Tend.row,
method = "linear")## Warning: adjust_rate: One or more of the timestamps for the rate(s) in 'x' do not lie between the timestamps for the 'by' and 'by2' background rates.
## Ensure this is correct. The adjustment value has been calculated regardless by extrapolating outside the background rates time window.## Warning: adjust_rate: background rates in 'by' and 'by2' differ in sign (i.e. one is +ve, one is -ve).
## Ensure this is correct. The 'linear' adjustment has been performed regardless.## adjust_rate: Rate adjustments applied using "linear" method.Converting units
apoly_cr.int_adj2 <- apoly_cr.int_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253) ## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.Plot curve
ggplot(as.data.frame(apoly_cr.int_adj2$summary), aes(x=row, y=rate.output*-1)) +
geom_point() +
stat_smooth(method = "lm", formula = y~poly(x, 2), color="red") +
theme_classic()
Rate filtering
apoly_rmr <- apoly_cr.int_adj2 |>
select_rate(method ="rsq", n=c(0.95,1)) |>
select_rate(method="lowest", n=6) |>
plot(type="full") |>
summary(export = TRUE)## select_rate: Selecting rates with rsq values between 0.95 and 1...## ----- Selection complete. 2 rate(s) removed, 19 rate(s) remaining -----## select_rate: Selecting lowest 6 *absolute* rate values...## ----- Selection complete. 13 rate(s) removed, 6 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: 14 1 298.8714 -0.01583195 0.979 NA 5211 5393 12669.04
## 2: 15 1 310.1747 -0.01599929 0.960 NA 5612 5794 13209.79
## 3: 17 1 317.6944 -0.01532933 0.961 NA 6412 6592 14290.21
## 4: 18 1 333.1337 -0.01578466 0.963 NA 6810 6954 14829.55
## 5: 20 1 316.5871 -0.01368371 0.950 NA 7530 7711 15910.30
## 6: 21 1 325.9582 -0.01379161 0.952 NA 7928 8108 16450.41
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 12914.46 97.762 94.282 -0.01583195 -0.002498924 -0.013333031 -0.013333031
## 2: 13455.42 98.253 94.480 -0.01599929 -0.002708259 -0.013291033 -0.013291033
## 3: 14535.37 98.190 94.512 -0.01532933 -0.003126338 -0.012202988 -0.012202988
## 4: 15074.95 98.445 94.724 -0.01578466 -0.003335133 -0.012449524 -0.012449524
## 5: 16154.94 98.238 95.118 -0.01368371 -0.003753283 -0.009930428 -0.009930428
## 6: 16695.55 98.566 95.150 -0.01379161 -0.003962426 -0.009829189 -0.009829189
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04565 0.0005137 NA 36 27 1.013253 -0.1422441
## 2: %Air sec 0.04565 0.0005137 NA 36 27 1.013253 -0.1417960
## 3: %Air sec 0.04565 0.0005137 NA 36 27 1.013253 -0.1301881
## 4: %Air sec 0.04565 0.0005137 NA 36 27 1.013253 -0.1328183
## 5: %Air sec 0.04565 0.0005137 NA 36 27 1.013253 -0.1059432
## 6: %Air sec 0.04565 0.0005137 NA 36 27 1.013253 -0.1048632
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -276.9010 NA mgO2/hr/kg -276.9010
## 2: -276.0288 NA mgO2/hr/kg -276.0288
## 3: -253.4322 NA mgO2/hr/kg -253.4322
## 4: -258.5523 NA mgO2/hr/kg -258.5523
## 5: -206.2356 NA mgO2/hr/kg -206.2356
## 6: -204.1331 NA mgO2/hr/kg -204.1331
## -----------------------------------------Results
results <- data.frame(Clutch = Clutch,
Replicate =Replicate,
Male=Male,
Female=Female,
Population = Population,
Tank = Tank,
Mass = mass,
Chamber = chamber,
System = system1,
Volume = chamber_vol,
Date_tested = Date_tested,
Date_analysed =Date_analysed,
Swim = Swim,
Salinity = salinity,
Temperature = as.numeric(unique(firesting2$temperature)),
Resting_kg = mean(apoly_rmr$rate.output*-1),
Resting = mean(apoly_rmr$rate.output*-1)*mass,
rsqrest =mean(apoly_rmr$rsq))
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 8 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.0005137 | ch1 | Asus | 0.04565 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 254.23 | 0.130598 | 0.9626 |
Maximum oxygen consumption
Data manipulation
Inspect file
## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 4.70
## -----------------------------------------
Subset data
cycle1.start <- Cycle_1.mmr[1,1]
cycle1.end <- tail(Cycle_1.mmr, n=1)[1,1]
cycle1.start.row <- which(firesting2_mmr$TIME == cycle1.start); cycle1.start## Warning in which(firesting2_mmr$TIME == cycle1.start): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="## Warning in which(firesting2_mmr$TIME == cycle1.end): Incompatible methods
## ("Ops.difftime", "Ops.data.frame") for "=="cycle1_data <- firesting2_mmr |>
subset_data(from = cycle1.start.row,
to = cycle1.end.row,
by = "row") ## subset_data: Multi-column dataset detected in input!
## subset_data is generally intended to subset data already passed through inspect(), or 2-column data frames where time and oxygen are in columns 1 and 2 respectively.
## Subsetting will proceed anyway based on this assumption, but please ensure you understand what you are doing.## inspect: Applying column default of 'time = 1'## inspect: Applying column default of 'oxygen = 2'## Warning: inspect: Time values are not evenly-spaced (numerically).## inspect: Data issues detected. For more information use print().##
## # print.inspect # -----------------------
## dTIME ch1
## numeric pass pass
## Inf/-Inf pass pass
## NA/NaN pass pass
## sequential pass -
## duplicated pass -
## evenly-spaced WARN -
##
## Uneven Time data locations (first 20 shown) in column: dTIME
## [1] 1 2 3 4 6 7 8 9 10 11 12 13 15 16 17 18 19 21 22 23
## Minimum and Maximum intervals in uneven Time data:
## [1] 1.32 1.41
## -----------------------------------------
Calculating MMR
## Warning: auto_rate: Multi-column dataset detected in input. Selecting first two columns by default.
## If these are not the intended data, inspect() or subset the data frame columns appropriately before running auto_rate()
##
## # summary.auto_rate # -------------------
##
## === Summary of Results by Highest Rate ===
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 348.3374 -0.04420749 0.9887977 NA 51 96 5702.99
## 2: NA 2 348.2615 -0.04419530 0.9887485 NA 50 95 5701.65
## 3: NA 3 347.5445 -0.04407080 0.9883762 NA 49 94 5700.32
## 4: NA 4 346.3223 -0.04385522 0.9873299 NA 52 97 5704.35
## 5: NA 5 345.0482 -0.04363231 0.9864485 NA 53 98 5705.71
## ---
## 174: NA 174 185.4308 -0.01623922 0.9708504 NA 155 200 5843.46
## 175: NA 175 185.1079 -0.01618459 0.9712072 NA 156 201 5844.79
## 176: NA 176 184.8326 -0.01613872 0.9716964 NA 159 204 5848.86
## 177: NA 177 184.8185 -0.01613616 0.9716706 NA 160 205 5850.21
## 178: NA 178 184.7947 -0.01613179 0.9718587 NA 157 202 5846.17
## endtime oxy endoxy rate
## 1: 5762.99 96.115 93.701 -0.04420749
## 2: 5761.65 96.201 93.672 -0.04419530
## 3: 5760.32 96.218 93.588 -0.04407080
## 4: 5764.35 96.165 93.707 -0.04385522
## 5: 5765.71 96.080 93.622 -0.04363231
## ---
## 174: 5903.46 90.692 89.558 -0.01623922
## 175: 5904.79 90.582 89.511 -0.01618459
## 176: 5908.86 90.514 89.473 -0.01613872
## 177: 5910.21 90.410 89.472 -0.01613616
## 178: 5906.17 90.579 89.464 -0.01613179
##
## Regressions : 178 | Results : 178 | Method : highest | Roll width : 60 | Roll type : time
## -----------------------------------------Adjusting
## adjust_rate: Rate adjustments applied using "mean" method.##
## # print.adjust_rate # -------------------
## NOTE: Consider the sign of the adjustment value when adjusting the rate.
##
## Adjustment was applied using the 'mean' method.
##
## Rank 1 of 178 adjusted rate(s):
## Rate : -0.04420749
## Adjustment : 0.0008872263
## Adjusted Rate : -0.04509472
##
## To see other results use 'pos' input.
## To see full results use summary().
## -----------------------------------------Converting units
mmr_adj2 <- mmr_adj |>
convert_rate(oxy.unit = "%Air",
time.unit = "secs",
output.unit = "mg/h/kg",
volume = chamber_vol,
mass = mass,
S = salinity,
t = as.numeric(unique(firesting2$temperature)),
P = 1.013253)## convert_rate: Object of class 'adjust_rate' detected. Converting all adjusted rates in '$rate.adjusted'.selecting rates
mmr_final <- mmr_adj2 |>
select_rate(method = "rsq", n=c(0.93,1)) |>
select_rate(method = "highest", n=1) |>
plot(type="full") |>
summary(export=TRUE)## select_rate: Selecting rates with rsq values between 0.93 and 1...## ----- Selection complete. 0 rate(s) removed, 178 rate(s) remaining -----## select_rate: Selecting highest 1 *absolute* rate values...## ----- Selection complete. 177 rate(s) removed, 1 rate(s) remaining -----##
## # plot.convert_rate # -------------------## plot.convert_rate: Plotting all rate(s)...
## -----------------------------------------
##
## # summary.convert_rate # ----------------
## Summary of all converted rates:
##
## rep rank intercept_b0 slope_b1 rsq density row endrow time
## 1: NA 1 348.3374 -0.04420749 0.9887977 NA 51 96 5702.99
## endtime oxy endoxy rate adjustment rate.adjusted rate.input
## 1: 5762.99 96.115 93.701 -0.04420749 0.0008872263 -0.04509472 -0.04509472
## oxy.unit time.unit volume mass area S t P rate.abs
## 1: %Air sec 0.04565 0.0005137 NA 36 27 1.013253 -0.4810951
## rate.m.spec rate.a.spec output.unit rate.output
## 1: -936.5294 NA mgO2/hr/kg -936.5294
## -----------------------------------------Results
results <- results |>
mutate(Max_kg = mmr_final$rate.output*-1,
Max = (mmr_final$rate.output*-1)*mass,
rsqmax =mmr_final$rsq,
AAS_kg = Max_kg - Resting_kg,
AAS = Max - Resting,
Notes=Notes,
True_resting="")
knitr::kable(results, "simple") | Clutch | Replicate | Male | Female | Population | Tank | Mass | Chamber | System | Volume | Date_tested | Date_analysed | Swim | Salinity | Temperature | Resting_kg | Resting | rsqrest | Max_kg | Max | rsqmax | AAS_kg | AAS | Notes | True_resting |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 48 | 8 | CSUD013 | CSUD017 | Sudbury reef | 300 | 0.0005137 | ch1 | Asus | 0.04565 | 2023-03-21 | 2024-06-27 | good/good | 36 | 27 | 254.23 | 0.130598 | 0.9626 | 936.5294 | 0.4810951 | 0.9887977 | 682.2994 | 0.3504972 | ||
| ### Expor | ting data |
## Rows: 333 Columns: 25
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (9): Male, Female, Population, Chamber, System, Date_tested, Swim, Note...
## dbl (16): Clutch, Replicate, Tank, Mass, Volume, Date_analysed, Salinity, Te...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.